TsAGI 


ANALYZING THE EFFECT OF SKIN POSTBUCKLING ON GENERAL 
STRESSES, STRAINS AND STABILITY OF COMPOSITE 

STRUCTURES 


Prepared by: 

G.N.Zamula 
K.M.Ierusalimsky 
G.S.Kalmykova 
V.P. Fomin 


1QQR 



2 


ABSTRACT 

The present paper is a final technical report within the NCCW- 1-233 research 
program (dated June 1, 1997) accomplished as a part of co-operation between United States' 
NASA and Russia's Goskomoboronprom in aeronautics, and continues similar NCCW-73 
and NCC- 1-233 programs accomplished in 1996 and 1997, respectively. 

The report concludes studies in two domains, "Analyzing the effect of skin 
postbuckling on general stresses and strains in a composite structure" and "Evaluating the 
effect of skin postbuckling behavior on general stability of a composite structure"; the work 
was fulfilled in compliance with NCC- 1-233 requirements (as of June 1, 1997). 

Also, the present studies may be regarded as a partial generalization of efforts in [1 , 2] 
conducted within the above programs in what concerns postbuckling behavior of composite 
structures. 
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INTRODUCTION 

The previous studies in [1, 2] paid major attention to postbuckling behavior of 
composite plates and built-up panels with consideration of nonuniform heating; the present 
effort deals mainly with the general strength and stability problem for a built-up structure in 
which some components have buckled. 

In thin-walled aircraft structures the skin is known to buckle locally under loads 
much less than ultimate loads. So there appears the nonlinear problem to determine general 
stresses and load-bearing capability of the structures with buckled skins (refer to [3]). For 
example, static strength of a wing box in bending is usually limited by general stability of the 
upper stiffened panel - in which the local buckling of skin cells between stiffeners may occur 
at a lower load than that of the panel. 

Upon skin buckling the panel is regarded as an anisotropic structure whose skin has 
some lowered (reduced) elastic characteristics that depend on the load level and are involved 
in the usual equations for computing orthotropic panel stiffnesses. Thus, two diverse 
problems have to be solved. 

At the first stage the load should be increased gradually and the subcritical general 
stresses of the structure with the skin buckled should be determined by utilizing a reliable 
procedure (for example, with a simplified beam-based model or the finite element methods.) 
During these operations the element stress resultants are nonlinear functions of the load, 
and their values are necessary in solving the general buckling problem at each value of the 
(varying) load. For built-up structures (including wings and fuselages) this second problem 
may be solved by analytical and/or numerical methods. Let us consider main aspects of 
these two problems taking into account features of composite structures. 

A report on the numerical method for evaluating stresses in a thin-walled structure 
with a buckled metal skin was delivered by the present authors in Stuttgart at the World 
Congress on Computational Mechanics, WCCM-II (see [3].) New reduction techniques and 
finite elements were used to transform the problem to classical iterations with variable 
elasticity parameters which are employed when calculating stresses and strains in 
nonlinearly elastic structures. Specific to the skin out of composites is the necessity to 
compute its secant reduced stiffnesses by utilizing some relatively simple relations - in 
particular, those derived in [1, 2], Thus, the analytical method of [3] applies to structures 
with composite elements. If required, the simplified solution to the problem on postbuckling 
behavior of a composite plate derived in [1, 2] may also be used to estimate the strength of 
each buckled skin cell with consideration of in-plane forces and variable bending 
deformation caused by out-of-plane deflection. This problem is discussed below in Part I. 

At each load value the second problem (on general instability of an isolated 
orthotropic panel or a cylindrical structure as a whole) may upon imposing some limitations 



be solved effectively by utilizing the numerical analysis method the authors developed in U 
However, the stability equations, unlike the stress/strain equations, are requiring t e 
tangential stiffnesses which relate increments in generalized resultants and strams. This 
forces analysts to compute the tangential stiffnesses (introduced in [1]) of the buckled 

composite skin. Relevant issues are considered below in Part II. 

In addition, both Parts pay much attention to algorithmic aspects and parametric 
analyses in which the features and the level of the influence of skin postbucklmg on global 
stresses and stability of composite and metal structures. Moreover, both Parts allow for 

mechanical and temperature-induced loads. 

It is clear that the proposed approach to allowance for the skin postbucklmg e ect 

on global stresses and stability of built-up structures is simplified and limited. However, the 

only real alternative thereto is a finite-element-based solution of a general geometrically 

nonlinear problem on very fine meshes, and this requires obtaining all local and general 

buckling shapes and bifurcation points for variable loads up to the structural failure loads. 

The latest finite-element software programs enable trying this way - but the result validity 

would always be in doubt. 
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Part I. Analyzing the effect of skin postbuckling on general stresses and 
strains in a composite structure 

1.1 Simplified (engineering) methods for evaluating nonlinear stresses and strains in a 
thin-walled structure with a buckled composite skin 

The influence of skin postbuckling on global stresses in a structure with a thin 
composite skin is very simple to consider by relying on engineering analysis methods that 
use a beam model suitable in dealing with long regular wing/fuselage sections. 

Let us briefly outline principles of the method. 

Consideration is given to a thin-walled section of a stiffened structure with a general 
cross-sectional shape; the global orthogonal coordinate system OXYZ is introduced so that the 

X and Y axes are in the cross-sectional plane and the Z axis is running along the section (see 
Fig. 1 . 1 .) In the cross-sectional plane the coordinate system origin is at an arbitrary point. 

Both the skin thickness 8 and the cross-sectional dimensions of longitudinal 

stiffeners are assumed to be small as compared with general cross-sectional dimensions of 
the structure. 

For the cross section contour (at Z = const ) we introduce the local coordinate 
system snz whose origin is at an arbitrary point, the s axis is tangential to the contour, the 
n axis is normal, and the z axis is parallel to the Z axis. The skin is made of a composite 
material with a symmetric layup of layer groups with various principal directions. 
Equivalent mechanical characteristics of the skin are determined on the basis of mechanical 
characteristics of orthotropic layers (E„, E 22 , G 12 , p 21 and p 12 = p 2 l E 22 /E n ), the 
directions of layer groups, and volumetric fractions of the layer groups. 

The equivalent mechanical characteristics E z , E s , = G sz and Msz correspond 
to the data for the orthotropic material. 

The analytical method relies on the following hypotheses: 

1. The cross-sectional contour (at z = const ) is not deformable in its plane, i.e., the relative 
contour deformation is zero (e s = 0). 

2. Strain in any cross section is described by the plane equation: 

£ z = E o + Gy* + 0 x y , ( 1 . 1 ) 

Herein, s 0 , 0y and 0 X are certain functions of z; the prime sign designates 
differentiation with respect to z. 

Now we use conventional stress/strain relations for an orthotropic skin 



to derive the normal stress a z profile over the sectional contour: 
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a * = r 


HszMzs 


( e o + 6yX + e x y) 


0 . 2 ) 


As usual, ct s = -p sz cr z . 

3. Since the skin thickness 5 is thin, the skin is assumed not to take bending, so its 

tangential stress x and normal stress cr z are distributed uniformly throughout the thickness 

direction. In this case it is convenient to introduce the normal stress resultant and the 

tangential stress flow by using, respectively, 

N 0 = o z 5 , q = x5 . (1.3) 

The flow q is tangential to the cross section contour. Similarly, stiffener sizes are small in 

relation to the cross section contour sizes; therefore, we are allowed to assume that the 

stiffener cross-sectional area f s is concentrated at the relevant contour point and that its 

force is P s = f s cr zs where is the stiffener stress. 

4. At each cross section z = const external loads are 


- bending moments M x (z), M y (z), 

-torsional moment M z (z), 

- axial force N z (z) , 

- shear forces Q x (z) and Q y (z), 

these are referred to the axes X , Y , and Z above. Positive directions for the external 


factors may be seen in Fig. 1.1. 

Note that these factors satisfy the relations, 


dM x 

dz 




(1.4) 


Evaluating the normal stresses 

Unknown values of s 0 , 0y and 0' x should be determined from equilibrium equations: 

n 

f a z 8ds + £ ct zs j f$ j - N z , 
i = 1 

% 

n 

fa z yds+ Zo 2si f si yj = M x ; ‘ (1.5) 

i = l 
n 

fa z xds+ Za zsi f si Xj = -M v , 
i = i 

Here, a z is the skin stress, 

cr zs j is the stress in an i-th stiffener, 

5 is the skin thickness, 

f si is the cross-sectional area of the i-th stiffener, and 

n is the total number of stiffeners in the section. 
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If stiffeners are rather numerous, the integrals can be estimated in a simple way by adding to 
each stiffener a half of the areas of the skin cells neighbouring with the stiffener (see Fig. 1.1.) 
In this case the integrals in ( 1 .5) become finite sums; using ( 1 .2) we have 

n 

E X F si e zi = N z , 
i=l 


ElF siyi e 21 = M y , 

i=l 

n 

ElF si x i£2i = -M y , 


( 1 . 6 ) 


i=I 


where 


E S i r si F S i + 2 j J^(i) r °J ^ ^ ^ s j ’ 


(1.7) 


Summation in (1.7) is for the set {co(i) } of panels neighbouring with the i-th node, 


r si = E S j/E, r 0 j = 




(i-p; 2j n s 2S j) E 

E is the Young's modulus of the material to which the characteristics of stiffeners and skin 
cells are reduced, 

E z j is the Young's modulus of a j-th cell, E|j = (pjE z j, 

<p] is the secant reduction coefficient for the j-th cell (with cp j < 1 if the cell buckles), 
u s and u s are secant Poisson's ratios for the j-th cell upon buckling, and 


SZ j 


ZS j 


bo; is the width of the j-th cell. 

Formulas for computing p s , p s . and <p * are presented in [1], 

S Z J ZS J 

As (1.1) postulates, e zi = e D +0’ y x i + G^yj, so we can re-write (1.6) as, 


E- 


F S x S y 

S I I 

^ x 1 x A xy 

S y l xy l y 


l 0 yj 


N 2 
\ = \ M, 


-M, 


(1.8) 


Here, F = lF s j is the generalized cross section area, • 
i=l 
n 

S x = 22 F s j y j is the generalized area moment for the x axis, 
i = l 
n 

S v = XF si Xj is the generalized area moment for the y axis, 
y i = i 

I x = ZF si yf is the generalized moment of inertia for the x axis, 
i = I 

n 0 

ly = Z F si x f is the generalized moment of inertia for the y axis, and 
y i=i 

n 

I xy = I yx = Z F si x , y | is the generalized product of inertia. 


i = I 
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The system (1.8) is nonlinear because the generalized sectional characteristics include the 
reduction factors <p] which depend on skin stresses (upon buckling) — and on unknown 

values of e 0 , 0' x and 0' y . 

Therefore, the normal stress in structural components should be determined by 
iterations with external loads being increased in increments. At a fixed external load (the 
vector (N z , M x , M y }) the analysis runs as follows. 

Initially, all reduction factors <pj are set to 1. Equation (1.7) is used to compute F si 

for all values of i and establish the matrix with generalized sectional characteristics. The 
system (1.8) is solved to evaluate e 0 , 0' x and 0' y . The formula (1.2) provides stresses in skin 

cells. If these stresses are greater than critical for at least one panel, the respective secant 

reduction factors cp] and secant Poisson's ratios are calculated by employing relations from [1]. 

The reduction factors and Poisson's ratio values are used to obtain F si from (1.7) and the 

matrix of generalized sectional characteristics. The system (1.8) is solved to again evaluate e 0 , 

0' and, 0' v . The successive iteration process does usually converge for 3 or 4 cycles, 
x y 


Determining shear stresses 

To determine the skin shear stress flow q we should consider the equilibrium 
condition for a skin element dsdz (Fig. 1.2) at N z = const: 


da 7 dq 

-^5 + ^ = 0 ‘ 
dz os 

Then the shear stress flow may be written as, 


q = 9 q + Qo( z )’ 


(1.9) 

( 1 . 10 ) 


do 


where q n = -I — ~ §ds and q 0 is an unknown closing flow. 
y - dz 


Now we employ the above discrete section model to write: 

q Q ,=-Z% L 5,s,=-EiF il |(s o+ 9;x, + e;y i ) = = -E|:F,,( £ ; + e; 5 <, + 0;v 1 ). HJi) 

Differentiate equations (1.8) with respect to z and involve (1.4) to obtain, 

( 1 . 12 ) 

Thus, for us to determine the derivatives e' 0 , 0 y and 0 X (present in (1.11)) we have 
equation (1.12) in which the matrix of generalized sectional characteristics is the same as 



F 

S x 

s y 



o ' 

E- 

S x 

lx 

1 xy 

K 


Qy 


S y 

^ xy 

] y 

k. 


[QxJ 
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that in (1.8), and the right-hand side vector contains shear forces. With this set of equations 
solved, the shear stress flow is calculated by (1.11) taking into consideration that 


F]=Z F si . , S yi =ZF sl 


Xi - 


(1.13) 


i=l i=l i=l 

are the generalized sectional area and the static moments of the cut part; these are computed 
allowing for sectional topology. The closing flow q G is established by satisfying the 

equilibrium equation for the torsional moment M z . 

If the sectional contour is simply connected, the flow q Q is described by the Bredt 


formula: 


in which the symbol Q denotes the doubled internal area. 

If, however, the section is multiply connected, the closing flows q Q j should be 

estimated by involving both the equilibrium equation and strain compatibility relations for 
edges of the contours; taken together, these provide the number of equations necessary for 
us to determine all flows q 0 j. 

In case the skin cells buckle under compression/shear, this circumstance is allowed 
for when estimating the reduction factors cpj by the method explained in [1]. Problems on 
normal stresses and shear flow are solved jointly at each iteration for determining reduction 


factors cp] and Poisson's ratios and Pszj with consideration of shear stress. 

This approach has been implemented in the composite structure section behavior 
analysis software. 

These means were utilized to investigate into the influence of postbuckling skin 
deformation on global stresses in a stiffened cylindrical structure. Two examples were treated 
in order to evaluate in numbers the significance of these factors in stress analyses. 

The first example is a circular cylinder with a 2540-mm diameter, loaded with a 
bending moment and a torsional moment, see Fig. 1.3. Distances b s between longitudinal 
stiffeners are identical throughout the contour, and the frame spacing £[ is 500 mm. The 
skin is a 2.64-mm thick composite stack comprising 22 layers whose stacking sequence is 
[0 5 /90 i/45 2 / -45 2 ] s . Equivalent skin data computed for the particular volume fractions 

and ply angles are as follows: E z =6682kg/sq. mm , E s = 3768 kg/sq. mm , 

=1434 kg/sq. mm, p.^0.3187. Similarly, skin ultimate stresses are o c [z =44.5kg/sq.mm 
in compression, a\ 7 = 31.2kg/sq.mm in tension, and x r2s =24.8 kg/sq. mm in shear. 



II 


Consideration is given to two versions of longitudinal stiffening (see Fig. 1.4); these 
differ in the stiffener spacing (b sl -lOOnrn. and b s2 =125tnnt) and the stiffener section 
area ( f sl = 200sq. mm and f s2 = 247sq. mm ). The sectional areas, f sl and n f s2 , are chosen 
such that the total stiffener areas for the two versions are equal. With this, the equ, valent 
stiffener data are as follows: E z =8346 kg/sq.mm, ultimate compressive stress 

55.5 kg/sq.mm, and ultimate tensile stress o'* - 39.0 kg/sq.mm . The structures are 
analyzed while proportionally increasing the external load vector: 

M x =tM x0 , M x =tM x0 , 

(where t is the increasing proportionality coefficient) for two load patterns: 


a) pure bending with M x0 = 5 • 1 0 6 kg • mm and M zo = 0 , 

b) bending and torsion with M x0 = 5 ■ 1 0 6 kg • mm and M zo = 2.55 • 1 0 6 kg • mm . 
Calculated stresses for these stiffening schemes and load conditions are shown in 

Figs 1.5 and 1.6 representing the lower stiffener stresses and mean stresses in the lowest 
(compressed) skin cell; the argument here is the load parameter t = M x /M xo . Figure 1.5 is 
to the load condition a); and Fig. 1 .6, to the load condition b). Let us first discuss the results 


for the structure loaded in bending (Fig. 1 .5.) 

Two straight lines (rays) which are running through the coordinate system origin and 

marked with circles demonstrate the linear analysis not allowing for skin reduction at the 
postbuckling stage. The upper ray reflects stiffener stresses, and the lower one, the skin stress. 
Solid fine lines represent stiffener stresses and the mean skin stress in case the skin reduction 
upon buckling is carried out. One can see that in the structure with the spacing b sl = 100 mm the 

lower cell buckles at t * 0.58 . With the load growing, the stiffener stresses increase, whereas 
mean skin stresses fall in comparison with linear analysis results. In the structure with the 
spacing b s2 =125 mm the skin gets buckled at an earlier time: t*0.38; in other respects the 

stiffener stresses and mean skin stresses behave like those characteristics of the former version. It 
should be mentioned that stresses notably re-distribute in comparison with the linear analysis 

and the difference increases as the postbuckling level gets higher. 

Moreover, the skin reduction amount varies in the vertical direction, so the neutral 

axis is shifted to tension-loaded panels and this circumstance, in turn, causes the compression- 
loaded structural part to take a greater force. Thus, the compression-loaded stiffeners are 
carrying additional forces because of not only the "early" skin buckling but also the neutral 

axis shift. 

The results for the case with bending and torsion are depicted in Fig. 1 .5. These are 
similar to those in Fig. 1 .6. Note that skin buckling stresses o z have diminished since the skin 


I r m Q 


f ulcn r^rrv in 


.nl^ne shpnr In thk 


rase the lower rvmel features the ratle - ft < 
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When assessing strength of a buckled skin cell, we must compute the total 
maximum stresses through the thickness which are comprising the middle-surface stress and 
the bending stress depending on local skin displacement w(x,y). These stresses may be 

computed by the method exposed in [1]. 

Figure 1 .7 depicts the dependence of the lower skin maximum stress (referred to the 
ultimate stress 5 ) on the parameter t = M x /M x0 for the case a). The total stress <r 2 at the 
point { s/b s =0.5,z = ^ x /2,h = 5/2} reach the breaking condition 
-at M x /M xo « 0.92 for the cell with b sl =100 mm and 

- at M x /M x0 * 0.98 for the cell with b s2 = 1 25 mm . 

Figure 1.8 shows a similar kind of dependence for the load case b). Torsion 
decreases the value of t = M x /M x0 at which the breaking stress is attained: 

- M x /M x0 *0.76 for b s2 = 125mm and - M x /M xo *0.82 for b sl = 100mm. 
Coordinates of the point with the hardest load depend on the external load. As for the 
above value of M x /M x0 , the breaking stress is attained at the point with coordinates 

s/b s * 0.3 ; z * 0.3^ x . 

The second example is a rectangular cylinder with the 2000-mm width and the 1000- 
mm height, loaded with bending or bending+torsion, see Fig. 1.9. The rib spacing £ { equals 
500 mm. Upper and lower panels and vertical side walls are 4.56-mm thick composite 
materials with 38 layers stacked as [0 10 / 90 3 /45 3 /-45 3 ] s - Equivalent skin data are as 

follows: E z =7383 kg/sq. mm , E s =3407 kg / sq. mm , 0^=131 1 kg/sq. mm , ^=0.3217. Skm 
ultimate stresses are o e h =49.2 kg/sq. mm in compression, a\ z = 34,5 kg/sq. mm in tension, 
and Tfo = 26.7 kg/sq. mm in shear. 

Consideration is given to two versions of longitudinal stiffening, see Fig. 1.10; these 
differ in the stiffener spacing (b s , =125 mm and b s2 =167 mm) and the stiffener section 
area ( f sl = 247 sq. mm and f s2 = 323 sq. mm ). The sectional areas are specified so that the 
total stiffener areas for both versions are equal. Equivalent stiffener data are the same as for 
the first example. 

Structures are analyzed for two load patterns: 

a) pure bending (with the upper panel being in compression and the lower panel subjected to 

tension) and 

b) bending + torsion 

while proportionally increasing the external load vector for M xo =3 10 kg- mm and 
M zo = 1.97 • 10 6 kg -mm. 
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Stresses calculated with and without allowance for nonlinearity are shown in Figs 
1.11 and 1.12; as in the case with the circular cylinder the skin and stiffener stresses are 
functions of the load parameter t. Figure 1.11 is for the pure bending, and Fig. 1.12, for 
bending + torsion. Figures 1.13 and 1.14 represent how the maximum stress in the compressed 
skin varies with t; here, the "circle" ( » ) sign denotes values o,/5 = 1 that correspond to loads 

under which the stress at the external skin surface become equal to the ultimate values. 

Qualitatively, the postbuckling deformation of the skin influences global stresses 
and strains in the same manner as in the circular structure. Note that the wing skin is 
usually thicker than the fuselage skin and is allowed to buckle under loads that exceed the 
in-service values. Therefore, the skin reduction influence on global stress/strain fields in the 
wing is less pronounced than that in the fuselage. In addition, wing skin reduction does not 
almost result in additional load for compression-loaded panels - because spar webs give a 
little contribution to the wing bending stiffness, and the forces taken by upper and lower 
panels are mainly defined by the bending moment and the section height. In the example 
above, the typical wing situation is modeled by the structure with b sl -125 mm and 

f sl = 247sq. mm , and the skin buckles at t « 0.93 . 


1.2. Nonlinear finite-element methods taking into account skin postbuckling behavior 

An overall algorithm for numerically analyzing strength and stability of a thin-walled 
structure with a buckled skin is depicted in Fig. 1.15. Here, t stands for any increasing 

parameter (it may be time) which is an argument for loads; 0 < cp S ^ 1 and 0 < cp l < 1 are 
secant and tangential reduction factors that evaluate degradation of elastic characteristics of 
the locally buckled skin; u is the value of the parameter at which the structure fails because of 

a) general buckling or 

b) destruction of primary components not buckled or 

c) destruction of the buckled (composite or metallic) skin. 

In the latter case we must pay attention not only to strength criteria with in-service loads but 

also to 

- requirements on allowable skin displacements and 

- the criterion of not reaching the yield strain (for metals) or of first-ply failure (for 
composites) under in-service loads. 

The problem is nonlinear because the skin reduction factors depend on skin 
stresses; note that, according to [1 , 2], the buckled skin is modeled by a stiffness-equivalent 
non-buckled orthotropic plate made out of a nonlinearly elastic material. 
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When analyzing global stresses of thin-walled structures by means of finite-element 
programs, the stiffened panels and walls are usually modeled with triangular and 
quadrilateral orthotropic plane-stress finite elements (Fig. 1.16a); nonlinearity is dealt with 

at each value of t by conducting additional iterations which refine <p s . Stiffeners 
surrounding a skin cell are usually spaced rather closely and deform in compatibility with it; 
therefore, the problem is below proposed to be solved by using new types of built-up finite 
elements comprising the skin and stiffeners; these are convenient in allowing for 
postbuckling behavior. In this case the composite skin orthotropy axes are assumed to 
coincide with stiffener directions; the skin and stiffeners may have different temperatures. 


Finite-elements 

Consider a thin-gage triangular or quadrilateral plate of thickness 5 orthogonally 
stiffened with ribs (in x- and y directions, at spacing , i 2 and cross-sectional areas F] , F 2 ). 
A skin material is assumed to be orthotropic, with principal axes running along the stiffeners. 
Local coordinate systems have x m , y m and x' m , y' m planes in a middle surface of the skin 
(Fig. 1.16b), the structure is located in a global coordinate system Oxyz as desired. A skin 
temperature field is considered uniform, of value T sk , whereas a stiffener temperature varies 
along the z - axis, the same law being assumed for all stiffeners of the same direction. To 
transform this model to an orthotropic material membrane, use can be made of the 
«orthotropic layer» concept (by Birger, Ref. 5) with «core coefficients)) K x , K y and K y for 

x' m and y' m — directions and for shear, respectively. Then, the elasticity relations for the 
layer may be written as 


e x =— 

E„ 


v y K y a y a y 

~ F +A x> 8 y = i 7 


x~-x~x . * T xyKx 


v,K,o, 


+ A y , y xy = 


GK, 


T yx K y 

GK„ 


0-15) 


or, in terms of stresses, 


° x l-v v vK K, 


-(£ x " A x ) + 


X y' 

E, 


v K E 


X y 


j 2 , _ . v xK-xEy 

CTy ~ l-v v v K V K„ y y) 1 - v v v„K v K 


(^X ^x)> 


(1.16) 


x y' 


GK 1 

Txy = ~K~ 


GK, 


’ Y xy * ^ yx 


x ' y “ x y 
’ ( V v M’ xv ’ ^ X ~ Mxv )■ 


T^r I XV’ ' y r* xy * x r* xy 


Here, A x = a x (T-T 0 ) , A y = a y (T-T 0 ) are the thermal expansions, T 0 is the initial 

temperature of panel. Integrating over periodical structural elements (Fig. 1.16c) and 
assuming the in-plane deformation, one can derive the following elasticity relations: 
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^x ^lx( e x £ xo)^"^2(^y ^ yo )’ 

Ny “f2^x — ^xo ) ^ly^y — ^yo)’ 
^xy "" ^yx — ^lyY xy» 


(1.17) 


where 


h — 6/2 


f.,= i 


E,K, 


h— 5/2 


-8/2 >- v x v y K x K y 


dz 


m 


L ly 


= J ^ 

*n 1-v.v 


h— 8/2 E V K K 

f _ | IL x v y Jv x iV y 


-6/2 ^ v x v y^x^y 


-dz 


m 


-6/2 1 "x T y^x J 
h— 6/2 

f, | GK r dz 


E„K y 

— dz , 

K V K„ m 


(1.18) 


- 8/2 


m 


are stiffnesses of the panel, 
flyTlx “ ^2 T ly 


e xo — 


My-fj 


e yo 


f|xjly "^ix 

Vi, - fi 


u u / 

r ix = I 


h " 5/2 E X K X (A X +v y K y A v ) 

-5/2 ' x v yK x K y 


h-5/2 


EvKv(A v "t" v x ^x^x) 


(1.19) 


y.- y“y/ „ f ^y^yv^y 1 y x— x-x 

dz m j T|y — J : 1 — ~ 


- 6/2 


1 v x v y^x^-y 


m 


1 v x v y^-x^-y 

are temperature strains and stress resultants. 

Eqs (1.15) to (1.19) are the same as formulae of Ref. 4, if the material used is 
isotropic, but they are convenient for FEM codes to use. Really, eqs (1.17) correspond to 
the typical matrix formulation: 

{N} = [D]({s} - {e 0 }) , (1.20) 


where 


[d] = 


f 2 

f .y 


0 

0 


0 0 f, 


ir- 


is the constitutive matrix, {N} — [N x N y N xy ] , {c} — [c x ®yYxyl * I^xo^yo^l are 

vectors of stress resultants, and initial strains, respectively. After dividing eq (1.20) by the 

skin thickness 5 and after comparing (1.20) to (1.16) at K x =K y = K y =1, one can see 

that, in place of the panel, there is a membrane of thickness 8, temperature T sk , average 

* 

stresses {cr} = [o x a y t xy ] T = 1/5{N) , made of an orthotropic material whose characteristics 
is principal directions ( x^ , y^ ) are E x , E y , G , v y , a x , a y , so that 

E. 


{5} = [D]({e} - {e 0 }), [D]=^[D] = 


E x V y 


l-v x v y l-v x v y 
E y v x E y 


l-v x v y 1-V x v y 
0 0 


0 

0 

G 


( 1 - 21 ) 


where {e„} = fa v (T„ fi -TJ a v (T nR -T,,) ol and 
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= - f| - 


E* = 


«* = 


p _£lx£ly_j2_ r _f - _ - h - E* 

r * r , G — f|y , Vy — , V v — „ — * v, ( — 

ii„ i] x t| x 


iy 

E xo - 

t _t ’ a >' 

*06 x o 


f, 


iy 


y E > 
y 


w f _ fix r _ *<y ? _ £ 2 ? _ £b_ 

T o6 -T n ’ lx 6 * ly 5 ’ 2 " 5 ’ lY " 6 ' 


(1.22) 


After integration in (1.18) and (1.19), formulae (1.22) can be used for specifying 
material properties in usual structural analyses with triangular and/or quadrilateral membrane 
FEs. However, it is important to note that the actual stress field is represented by eqs (1.16) 
instead of (1.21). In particular, a skin experiences, as rule, additional local stress due to 
restraint from stiffeners of lower temperature; this circumstance become important while 
analyzing whether the skin can buckle (and at reducing the skin, if any). Hereinafter, let’s use a 
simplifying assumption of a small values of K y «1, and v y K y «1, v x K x «1 for all 


parts of the panel under study. Such an assumption is usually valid for typical airframes. 
Then, the skin is characterized by K x = K y = K y = 1 . Eqs (1.18) will be transformed: 


f - 

r lx - 


E y 5 E 1 F 1 


+ -b L - f 'y=7TT7Tr + -T i . f 2=JT^-. f i, = G«> 0 23) 


E V S E?F, 


E v v v 5 


l-V x Vy £, 


l-v Y v v £ 2 


X Y y ^ l * r X Y y 

where E x , E y , G , v y are the elastic constants of skin at temperature T sk , 


E ,=pJ E ,dF, E 2 |E,dF 


= 1 


are the averaged moduli of stiffeners. By parity of reasoning with eq. (1.19), obtain 
E v 5 E|F, _ _ 1 

*\ 

E.,5 so . n 0-24) 


T i* = + v y“y)( T o6 -T 0 ) + -)-La„ A, = Je^CT, -TJdF, 


\ r \ 


T, v = 


1 — V V 

1 V x Vy 


E F 1 

(a y + v x a x )(T o6 -T 0 ) + -j-^A 2 , A 2 = =^r jE 2 a 2 (T 2 -T 0 )dF, 


where Aj and A 2 are the averaged thermal elongations of stiffeners. Eqs ( 1 . 1 6) give stresses 


E x (l-v y v x )_ E (v - v ) 


<*x = 


a x + — 


CT, ,-T 


E x (l-v x v y ) * E y (l - v x v y ) y l-v x v y 


^ xo yo 


a = =- 


E y(l ~ V x V y ) _ E y (v x -V x )_ 


■CT„+ ~ 


a -7 


y E y (1 - v x v y ) ” y E x (l-v x v y )^ 1- v x v y 


-[(a x +v y a y )(T o6 -T 0 ) 

[(a y + v x a x )(T O 6-T 0 )-e y0 - v x e, 


,(1.25) 


G 
G 

for a skin, 


T = ■=■ T = T 
u xy q w xy 1 xy 


Ej Ej Vy p 

" £ a x " £ “ E||a|(Tj — T 0 ) — e X o]» T xy = 

c y - g CT y “ £ °x-E 2 [a 2 (T 2 -T o) - eyo ], t yx = 0 


(1.26) 


for stiffeners. 
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It can be seen from ( 1 23)-(1 .26) that both determination of skin stresses and 
reduction of skin require only a few data: skin parameters and integral charactenst.cs of 

stiffeners, namely, F, , (\, E,, A,,F 2 , t-i, E 2 , A 2 . 

To perform a strength analysis of such finite elements, one should apply usual 
relations or the Finite Element Method, Ref, 5; in addition, an elastic and thermal 
performance of the skin and the stiffeners can be specified. Hereinafter, the parameters 
involved are assumed to be the same over all skin and over all stiffeners of the same direct, on. 

A .rtanunlar stiffened panel (Fig. 1 . 1 6b) carries in-plane axial loads and shear. The 
displacement field in the element is assumed to be linear. For the plane-stress element with a 
uniform thickness 8 the stiffness matrix is, 

[K] = 5 J[B] T [D][B]dA , (127) 

A 

Oere 

[D] is the matrix with reduced element stiffnesses referred to the local { x m , y m } 
coordinate system and derived from [D] (see (1.21)) by transforming the coordinates; 

[B] is the usual constant matrix for computing the strains {s} = [B] {q} . 

Upon integrating the relation ( 1 .27) the explicit formulas for the element stiffnesses 

referred to the local coordinate system appear. 

Equilibrium equations taking into account initial thermal strains are, 

[K]{q} = {F} + {F} £o 

Here 

{F} is the nodal force vector due to the external loads and 

|F> is the initial thermal force vector which may be evaluated as 

l t E 0 

{F} Eo =5 J[B] T [D]{e 0 }dA 

The FEM problem is then solved to determine nodal displacements {q}, strains {e} , 

thereafter, (1.21) is utilized to calculate average stresses, and relations (1.25) and (1.26) are 

used to compute the stress field for the panel. 

A quadrilateral stiffened panel (Fig. 1.16b) is to be modelled by a simple finite element 

whose thickness is that of skin, and thermo-mechanical characteristics (1.23), (1.24) take 
account of performance of both skin and stringers. A stiffness matrix can be derived as 
suggested in Ref. [6], That is, the quadrilateral structure should be separated into four triangles, 
and displacements of the additional internal node should be eliminated. Each triangle has a 
linear displacement field and a uniform stress state. For each triangle, a structural analysis 
system can compute a stiffness matrix, strains, and stresses by means of the above formulae. 

Unfortunately, such FEs produce significant errors if used in areas of severe stress 

ararlients. sav. in models of spar walls subjected to bendine. 
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An accuracy can be improved through introducing the quadratic incompatible 
shape functions, their coefficients being found under conditions of minimum strain energy 
of the finite element. Ref 7. A new FE — quadrilateral stiffened wait - has been a good 
complement for the existing FEs of our structural analysis system, ft represents a futher 
improvement of the FE proposed in Ref. 7 and allows for additional internal structures, 
orthotropic material, and temperature gradients. The stiffness matrix components can be 
computed using the Gauss-Legendre four-point integration formulae, Reft 6; eq. (1.21) is a 
source of the matrix of elastic constants. The finite elements proposed can be developed 
further to include different structural features. 


Solving the nonlinear problem 

As stated above, consideration of plastic deformation and buckling of the skm 
makes us regard stiffnesses of the skin (and the panel as a whole) as functions of an 
instantaneous level of stresses and strains: 

[DMD«e},(a)», fe 0 } = {e 0 ({s),{o})), (L28) 

where la! = [o,ct v t„] t is the skin stress vector in accordance with (1.25). 

The nonlinear problem may be solved on the basis of the FEM and the proposed 

finite elements by resorting to iterations. At each stage we should analyze linearly the 
stresses and strains by using the stiffness matrices and temperature forces which result from 
strains, stresses and secant reduction factors obtained at a previous stage. In this case the 

reduced stiffnesses of a buckled composite skin (E^, E’ y , G s , v*, v s y = v s x E'/E* ) with 
and without temperature variation consideration depend on the "supercritical strain" of a 
particular finite element, and one should rely on relations of [2] and [1), respectively. 

At the initial iteration we solve the usual elastic problem for the structure with no 
skin buckling. Iterations stop when difference of stresses as obtained at two sequential 
iterations becomes rather small. In particular, if an unbuckled metal skin is involved m 
plastic deformation the present iterative process coincides with the traditional variable 
stiffness method [6] for physically nonlinear problems. The new method has been 
implemented in a computer program and validated by using example problems for which 
theoretical and numerical solutions have been provided elsewhere. New numerical analysis 

results are rather verisimilar and close to experimental data. 

One knows two methods for finite-element modeling of a stiffened structure with a 

buckled skin. The first method is to prepare a rather dense mesh whose characteristic step 
sizes and i 2 are close to the stiffener spacing, to use the rod elements for modeling of 
stiffeners, and to introduce unstiffened plane-stress plates (with F, = F 2 =0 in (1.23) and 
( \ ^411 for modeling of skin cells. The second method is suitable for analyzing structures 
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with numerous stiffeners; here, use is made of the built-up finite elements described above, 
whose dimensions are much greater than and f 2 ; this enables analysts to increase the 
characteristic element dimensions to the distance between primary load-bearing components 
such as spars and frames. Of course, the second approach requires much less computation 
time and computer memory space and would be preferred. 

Consider an example: nonlinear analysis of an edge-clamped cylindrical stiffened 
shell loaded with a transverse shear force. Structural dimensions and the value of the force 
may be seen in Fig. 1.17. This example is typical of real fuselages. Results of experiments 
with the present shell are reported in [8], The skin is very thin (8 =0.51 mm), and its cells are 
subjected to various combinations of normal and shear stresses which lead the skin to 
buckling nonuniform over the cross section contour. 

Analyses were carried out by using the nonlinear FEM program above. Introduced 
as nodal points of the finite-element mesh were points of intersection of longitudinal and 
transverse stiffeners. Stringers and frames were modeled with beam-type finite elements, and 
the skin, with flat quadrilateral membranes. Thus, the finite element for the skin models a 
real skin cell between neighboring stringers and frames. The FE model includes 126 nodal 
points with the total of 630 unknown displacements. At each iteration the program 
corrected the stiffnesses of element that model the shell skin: initial elastic properties of the 
skin material were multiplied by current values of reduction factors. The reduction factors 
for the skin under multiaxial load were computed by employing the method of [1], 

Figure 1.18 depicts stress diagrams for the skin and stiffeners at the clamped end. 
Also, one can see experimentally obtained results from [8]. Comparison for the compression- 
loaded domain shows that data of experiments and the nonlinear analysis are in a rather good 
agreement. 

Figure 1.19 represents diagrams of normal stresses as computed for stiffeners at the 
clamped end by using linear and nonlinear theories. It is clear that stresses differ much, 
especially in the compression-loaded domain. As compared with the linear calculation, the 
mean normal stress in the skin over the severest compressjon area has decreased by a factor 
of over 2.3, whereas stiffener stresses in the same area have been higher by 23%. 
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Part 2. Evaluating the effect of skin postbuckling behavior on general 
stability of a composite structure 


2.1. Tangential stiffnesses of stiffened structure upon general 
buckling 


In Part 1, we dealt with analysis of stresses and strains in a stiffened structure whose 
thin composite skin carries loads while being buckled. At a certain value of the external load 
global instability may occur. At that moment both the stresses and strains show small 
increments which may be regarded as variations to initial (subcritical) values for the 
structure with nonlinear dependence of strains on skin stresses. 

In [1] the nonlinear dependence for a composite plate has been written as, 

Py Pv Py Px 9^. x v 

e x =^--p yx -^ + B e (p,s), e y = — -p xy — + a a B e (p,s), 

b x b y t y fcx (2.1) 

y = 7 + B 0 (p,s), 

°xy 


where a a = a / b and p = p x + a 2 p y . 

The relations (2.1) were obtained by assuming proportional loading — that is, at 
p = p* • t and s = s* • t , where t is the load parameter (t>l ) and { p*, s* } is the critical load 
combination (refer to Fig. 2.1). Functions B e (p,s) and B 0 (p,s) in (2.1) represent the fact 

that upon local buckling the initially orthotropic plate becomes anisotropic at s * 0 . 
Variations to the initial relations (2.1) look like this: 

5e x = a[|5p x + a| 2 8p y + a[ 3 8s, 

8e y =a 21 8p x +a 22 8p y +a 23 8s, (2.2) 

5y =a 31 5p x + a 32 8p y +a^ 3 8s, 


Here, 


a t L + gBe 

"~E X + dp 
t 1 4 <5B e 

a22= i: +aa “^’ 


a 12 = a 2l ~ 


Pxy 2 

— — + a 


3p 5 


a‘ 

3,3 ■ ds ’ 


a 33 - 


1 dB 
+ 


e 


a 3l = 


7 

dBe_ 
dp ’ 


G xv 


ds ’ 


t 2 ^B e 

a 23=« 


ds ’ 


a 32 =a 


-> 5B f 


dp ’ 


Potentiality conditions require that a [3 = a^j and a 23 = a 32 i th' s means that 

dB c _ dB 0 
ds dp 
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The latter condition will be met if the loading is proportional and the ratio of 
supercritical out-of-plane displacements, vp = f 22 /ft i » d Qes n0 * depend on the 
"supercriticality index" t; i.e., 

_ * 
s 


V = 


4-p 


(2.3) 


where p* = p / p x0 and s = s /s 0 . In this case 


and 


2 * 

B e = — y f] 2 ] (1 + 4y 2 ) , B 0 = 4^-*r B e vp , 
8a Sq 


ci p - 1 + 4\ps P - J_ 

In— i 

1 _ p x0 s 0 


(2.4) 




Expressions for C 7 and C 8 are provided in [1]. 

Let us symbolize derivatives of B e and B 0 with respect to p and s as follows: 


3B„ 


3B P 


dB f 


Bg p " ap C ’ Be s ds ’ Be p dp ’ 


d B e 
B fl , = — - 
e - s as 


(2.5) 


Then we take into account (2.3) and (2.4) to obtain: 


1 n 2 


Be -p _ Co 2 1 


1 -t-4vp J 


Pxo 8a^ _^c 7 -h/ z C 8 


1 n 

’> B e,s — * 


1+4^ 


R _ J_ PxO_ 

“e.s “ * * 2 

So s 0 a 


2 


1 + 4vp " 


■v 2 . 


So 2a^ ,i ( 2p 


•V = B e, P 


(2.6) 
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fC; 


With (2.4) and (2.5) written and the stability boundary equation being 
p* +4\|/s* =1, the original expressions (2.1) may be re-written in the form usual in studies 
on anisotropic materials: 

e x = a ?lPx + a l2Py + a i3 s > 
e y ~ a 2 1 P x + a 22Py + a 23 S ’ 

Y =a 3lPx +a 32Py + a 33 s > 


(2.7) 


E v + Bep 


a,, = 




1- 


, a j — a 2 1 ~ — “t - & o B 


Pxy 

E v 


2 

a ^e,p 


f + \ 

i-— 

k p j 



( 



f 

= B e.s 

.-i- 

> a 23 ~ a 32 

“ a a ®e,s 

S 

1- — 


l S ) 



l s J 


a 33 “ 


'xy V 




. S 


1- — 

S J 

\ 




B 0,s B c.s * 


e,s 


S J 


where 
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However, analyses are easier to perform upon transforming the relation (2.1) to the 
version corresponding to a quasi-isotropic materia, with variable elasticity parameters: 


e ’nr + 


i . iO 

x P t 


e v =1 - 


Px + 


Pxy 2 

J. /V 


Pyx 2 

: r* 


+ a 


Py 


+ a 


Px + 


+ a 


4 2e 

P ) 


( 2 . 8 ) 


1 Be 

Y= l^T s 


s. 


. .. c,.,.., «,* (C , 5 ) and the secant Poisson's ratio, the formulas 

As for secant reduction factors (<Px > y » ^y ' 

derived in [1] are valid. . 

Also, (2.2) may be re-written to the form usual in studies on orthotropic materials 

with variable elasticity parameters. Loading is assumed to be proportional, so 8p x = P x 5‘ • 

5p = p' St , and 5s - s‘ St ; for p* * 0 and s' * 0 we obtain: 


* \ 


* X 


8e x = a i , +a} 3 


5p x + 1 2 + a a a i3 p * y 


P ) 


6 py , 


5e y = 


s 

P 


a| 2 + a a a 13 ~y 5p x + l a 22 + a a a 13 


4 an — r Up, 

P ) 


(2-9) 


6y = 


^ a 33 + a{ 3 j 5s- 


Then the stress increment expressions which are inverse to (2.9) become 

+ + ^ '* = G " 5r ' 


1 Pxy P 


xyr-yx 


where 


a} 2 + <*a a i3 * 
p - 


E l = 

J-'x 


T , Pyx = 


t t 0 

a n + a 13 * 
P 


a}i + a i3 * 

P 2 


E l = 

E.y 


G l xy = 


t 2 „ t _ 

a i2 + a a a i3 * 

P 


t 4„t _L_ 
a 2 ? + a a a 13 , 

X " P - 


1 


t 4 t a 
a? 2 + a a a 13 * 

P 


t _t P_ 

a 33 + a i3 * 
s 


(2.10) 


We proceed from relations E x =E x cp x , E l y -E y (p y , and G xy Gxy^y’ releva 
reduction factors may be determined from (2.10). 
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For increments of stress resultants and moments upon global buckling the 
following equations are valid: 

N x = |8p x dz, N y = j5p y dz, N xy = J8sdz, 


M x = j8p x zdz, M y = j5p y zdz, M xy = J5szdz. 

h h h 


Now we can employ the Kirchhoff-Love hypothesis and the Shanley hypothesis on 
loading continuation at bifurcation; with this, the latter equations provide the following 
expressions for reduced tangential stiffnesses of the skin: 


Bfio = 
d!io = 


E‘h 


1 - FxX-x 

E l x h 3 


B 220 - 


E y h 


1 -FxvF 


B t _ G t . B t _ ^U E x h 

t t ’ B 330 -°xy n > ”l20 — 


;> D220 - 


xyf^yx 

E y h 3 


I FxyFyx 


12(l ~ Fxy^yx)* ^ 12(l VxyHyx)’ 


D33O “ 


G‘h 3 

, dU-m* d|, 0 . 


( 2 . 11 ) 
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If the skin is stiffened regularly with longitudinal and transverse elements (having 
spacings b s and b f , respectively) then the tangential stiffnesses of the stiffened panel are 
similar to (2.1 1): 


Bf, =E S 


rtsH+ ^. 




» B 22 = E f r tf h + ^J, B 33 -(^xyh +Gsfh sf j, 

B f 2 =B f 20 > B 3 = B330 + Bf 20 > 


Df, =E C 


L ts 


r h_ 3 

12 


+ hh 


+ ^(zs-h,) 2 



r , 

fh 3 -,1 

D22 = E f 

1 

Ti + hh2 

V lz J 


b s b s 




( 2 . 12 ) 


D33 = G 


t 

xy 6 


, , . G S I S P Grlf 

D\=2D^ 3 o+DU+ s s ■ r 


'120 


2b . 


+ 


2b 


f 


where 


El 


E* 


Trc = 


E s (l FxyFyx) Ef |l PxyFyx) 


E S ,E|-,G S , and Gf are elastic moduli and shear moduli of the longitudinal and 
transverse elements, 

P P 

F s , Ff , I s , If , I s , and l f are cross-sectional areas, central moments of inertia, and 
the product of inertia of the longitudinal and transverse elements, respectively, 
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z s and z r are distances from the skin midsurface to centroids of the longitudinal 
and transverse elements, respectively, 

h, and h 2 are distances from the skin midsurface to centroids of the longitudinal 

and transverse elements with relevant skin strips, respectively, and 

G sf • h sf is the shear stiffness of the frames/spars with no skin; this value depends 

on a method for joining the longitudinal and transverse frames/spars. 

These equations for secant and tangential stiffness parameters are used hereinafter 
to analyze subcritical stresses and strains and general stability of the stiffened structure. In 
the case with no skin buckling the equations become the traditional expressions for 

stiffnesses of a stiffened panel. 


2.2. Stability of regularly stiffened cylindrical and flat composite panels and shells 
in the case of uniform loading with compression and shear 

Consideration now is given to a stiffened cylindrical panel whose radius is R and 

planforra dimensions are A and B (see Fig. 2.2). The panel is regularly stiffened with 

longitudinal and transverse elements. The cross-sectional area and the spacing of 

longitudinal element are F s and b s , respectively; and those of transverse elements are F, 

and b r , respectively. It is assumed that the panel skin is a thin composite with layup 

sequence being symmetric with respect to the skin midsurface. Reduced mechanical 

characteristics of the skin correspond to characteristics of a homogeneous orthotropic 

r; t: n n (F u = E u ) The skin thickness is h. Dimensions 

material— E x , E y , G xy , p xy , P yx -^y^xy^- 1,16 

b s and b f are much lower than A and B, so the skin cell may be assumed to be flat. 

Let us use u, v, and w to denote projections (on axes x, y, and z, respectively) of a 
small displacement of a panel point. 

For a shallow panel either stiffened or unstiffened being uniformly loaded with 
stress resultants N° x , N° y , and N° Iy (with no moments applied) and unbuckled, we can write 
the following linear equations: 

[l]{u} = 0, (2J3) 

Here, [I] is a symmetric 3 by 3 matrix whose components are differential operators 
Ly for stability equations (i,j = 1,2, 3) and 
{U} is the displacement vector with components u, v, and w. 
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The differential operators in (2. 1 3) have the following form (as provided in [9]): 


Ln=Bf,^ + BL^, L n -L 2 l -B\ dx ^. 


L 1 3 — L^I - — 


31 


R dx 


, a 2 . d 2 

» ^22 = B 22^7 + B 33 . 2 * 


L 33 = D| 


11 dx 4 


+ 2D\ 


3 dx 2 dy 2 
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+ D22 
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^d 4 


2 d L 1 
as 4 R 2 ds 2 R 4 


L 23 - L 32 - 
B 22 


Bkl 

R ds’ 


(2.14) 


R 


-N 


°— - 2 N° 

X - 2 Y 

dx 


-N 


0 

y dy 2 


xy dxdy 

In the case of a panel whose edges x = 0 , x = A , y = 0 , and y = B are freely 
supported, the solution to (2.13) may be written as 


00 00 


nrmx . nrcy 

u=X lA mn cos— — sin 


B ’ 


n=l m=l 

® ® mrcx . nrty 

w=I ZC mn sin — sin 


® ® mnx nny 

v=I lB mn sin— -cos 

n=l m=l A 


B ’ 


(2.15) 


B 


n=l m=l A 

Upon substituting relations (2.15) into (2.13) and resorting to the Bubnov-Galerkin 
method we obtain an infinite homogeneous algebraic system of linear equations for the 
constants C r 


'mn * 


( c 11 c 23 ~ C 1 3 C 22 ) c 3 1 + ( c 13 c 21 ~ C 1 l c 23 ) c 32 + 


c 1 1 c 22 _c 12 


C 33 ~Pg 


(2.16) 


32amns 


r s 


n 


H-ZI C n — 2 , ; 2 . 

2 r S (m 2 - r 2 )(n 2 -s ) 


where m, n, r, s = 1,2, 3, ...00 and the values of m+r and n+s are odd numbers, 
c n =m 2 +b, 2 n 2 a 2 , c 12 = c 2 i = b 3 mna , c, 3 =c 31 =b )2 mp, 


c 23 = c 32 = b 22 naP , c 33 = 7i J 


dnm 4 +2d 3 m 2 n 2 a 2 +d 22 (n 2 a 2 -p 2 ) 


+ b 22 p , 


(2.17) 


•vrO isj 1 
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m + — ^n a 
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tO 
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Bf, 


N 

1 Y xy 
S p t > 

8 Bf, 


_ A _A_ h . P ‘j 

P ~-«’ ,J Bf, ’ 8 Bf,A 2 ' 


B 


The stress resultants N$ , Nj , and Nj y are assumed to be proportional to a single 
parameter t; one could limit the system (2.16) with a certain set of equations (it suffices to 
adopt m and n of 6 to 8); with this, it is easy to establish the critical value t = t which zeros 
the system (2. 1 6) determinant. 
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The global stability analysis is performed in compliance with the algorithm depicted 
in Fig. 1.15: 

the load parameter is incremented from zero in steps of At , 

_ at each load step the iterative procedure is run in order to determine subcritical 

stresses and strains in both the skin and stiffeners, 

— if the skin is to buckle at the current load level, then relations (2.10) and (2.12) 
are employed to establish tangential stiffnesses of the panel and their respective determinant 

of the system (2.16); 

— if the system (2.16) determinant changes its sign at a certain load level, then the 

critical value of t' should be refined to the accuracy required. 

The problem of general stability of a stiffened circular cylindrical shell (Fig. 2.3) is 

solved in a similar way. A particular solution to the system (2. 1 3) is represented by, 

u = A cos(X m x - nip) , v = Bcos(x„x-n<p), w = Csin(3. m x - n<p), (2.18) 

where x = x / 1 , cp = s / R , 

i is a shell half-length, 

5 is a cross section point coordinate measured along the circle, 

X m and n are wave parameters. 

By substituting (2.18) into (2.13) at a fixed value of n and performing the necessary 
mathematical manipulations we arrive at the characteristic algebraic equation of the eighth 
order with respect to X m ; roots of this equation taken together with the eight constants 
C (m = 1, 2, . . . , 8) make it possible to write the general solution to (2. 1 3). 

w = ZC m sin(^ m x-ncp). (219) 

m=l 

Constants A m and B m are functions of C m . 

Boundary conditions for shell ends { x = ±( } are the following. 

2 

w _ = N x = v = 0 in the case of free support at the ends, 

dx 2 

w _ — _ N = v = 0 in the case of clamping of the ends. 
dx x 

These allow one to formulate 8 linear homogeneous algebraic equations for 

constants C m . Coefficients of these equations depend on shell end loads N x , N s , and N xs . 
The condition of existence of a nontrivial solution to the system provides the critical load 

parameter value t* and the critical load (N x , N s , N xs } of global buckling. 

Let us delve into particular problems with stiffened structures whose skins buckle 
and carry loads upon buckling until the structure undergoes global buckling. 
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The initiatory problem is on a flat (1/R=0) stiffened panel with planform 
dimensions A=480 mm and B=400 mm. The skin with the thickness of 2.64 mm is made out 
of a composite with the following equivalent characteristics: E x = 6682 kg/sq. mm , 

E y = 3768 kg/sq. mm , G xy = 1434 kg/sq. mm , p xy =0.32. Ultimate stresses of the skin are 

a x = 44.5 kg / sq. mm in compression, a x = 31.2 kg/sq.mm in tension, and 

x sh =24.8 kg /sq.mm in shear. Longitudinal and transverse stiffeners are rectangular 

shapes with cross-sectional sizes 20 by 5 mm and the spacing steps b s =100 mm and 
b f = 120mm . Stiffener material characteristics are as follows: the Young modulus 

£=8346 kg/sq.mm, the compressive strength o c = 55.5 kg/sq. mm , and the tensile strength 


a 1 = 39.0 kg/sq. mm . 

Critical stresses of local buckling of a skin cell between stiffeners in the case of 
uniaxial loading are p* x0 = 10.87 kg/sq.mm , p* y0 =7.55 kg/sq. mm and sj=2 1.29 kg/sq.mm . 

Figures 2.4 and 2.5 show the analytical results. More particularly, Fig. 2.4represents 
the global buckling boundary curve in the plane {N x , N xy } (note that N y =0); it defines 

the critical load combination { N x , N xy } (shown by the solid line) and has been drawn with 
consideration of skin postbuckling behavior (until the panel becomes buckled) for the 
following load combinations: 


N 


xy 


N 


= 0, 0.2, 0.4, 0.6, 0.8, 1.0 and — ^- = 1.0, 0.8, 0.6, 0.4, 0.2, 0. 


N x N xy 


Also, one sees here a similar kind of boundary curves computed for the same panel 


without assuming the skin to carry loads upon local buckling (the dashed line). By 
comparing these curves it is clear that the global stability analysis with no correction for 
load-carrying capability of the buckled skin may be misleading in its overestimating the 
global buckling load by 25% to 60%. 

Figure 2.5 represents stresses in components of the stiffened panel being in 
subcritical state, immediately before global buckling; consideration is given to all the above- 
mentioned combinations of shear and compressive loads at the edges. Here, solid lines are 
for component stresses if the analysis takes into account load-carrying capability of the 
buckled skin, and dashed lines, for component stresses calculated with no account of load- 


carrying capability of the buckled skin. In case the shear load is relatively low the most 
notable difference is between mean stresses in the buckled skin, p x , and the nominal skin 
stress a x if the analysis does not allow for early skin buckling. When the shear load is great, 
one could reveal additional compression of both longitudinal and transverse stiffeners which 
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is caused by appearance of oblique waves in the skin see ct s and Of in Fig. 2.5 for 
N? -> 0 at Nj y * 0 . 

However, it should be noticed that, when the shear load is increased, the mean shear 
stress s_ in the skin grows rather rapidly and reachesthe ultimate level (x « 25 kg / sq. mm) 

at N xy =72 kg /mm . Thus, if the shear stress resultant N xy is higher than 72 kg/mm, we 

should expect the skin to break at loads less than the global instability load. 

Consider also analysis of global stability of a stiffened cylindrical panel with radius 
R=1270 mm. Panel planform dimensions, longitudinal and transverse stiffener shape and 
sizes, skin thickness and material characteristics are the same as those of the flat plate above. 

Figure 2.4 demonstrates the panel buckling boundary computed by the present 
method. Clearly, the panel curvature does notably improve general stability. Calculation for 
the same panel with no allowance for buckled skin stiffness provides results which are 

greatly overstated — by 20% to 40%. 

Consider now the problem of stability of stiffened circular cylindrical shell subjected 
to torsion. The shell has the following main parameters: curvature radius R=750 mm, overall 
length 2^ = 21 00 mm, skin thickness h=l mm. The skin is regularly stiffened with stringers 

(with the spacing b s =118 mm) and frames (with the spacing b f =150 mm). The stringers 
are rectangular 8 mm by 20 mm shapes, and frames are angles with leg sizes of 15 mm and 
35 mm and the web thickness of 2.0 mm. The material of both the skin and stiffeners is 
isotropic, E=7200 kg/sq.mm and p = 0.33 . 

Three identical metal shells have been manufactured and tested under torsion. At 
skin stresses x * 4.8 kg/sq. mm the skin had buckled and remained in this state until the shell 
failed due to its global buckling. Figure 2.6 demonstrates one of the shells in a subcritical 
state (as for global buckling). One can readily see local waves in skin cells. Here, the 
"supercriticality level" is approximately 3. Figure 2.7 represents global buckling of the shell; 
spiral waves involve stringers and frames into large deformation. The experiments provided 
the following values of critical shear for the three shells tested: 

Nxy(i) =14.2 kg/mm, N xy(2 ) =13.8 kg/mm, N xy(3) =13.7 kg/mm. 

Our analyses with postbuckling skin stiffness reduction suggest N xy =15 kg/mm . 

Thus, the experimental data are at the level of 91% - 95% of the theory. In contrast, analysis 
of the same shells with no allowance for stiffness reduction overstates the critical load by 

approximately 35%. 
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2.3 Analyzing the stability of a built-up cylindrical structure with consideration 
of global and local buckling, nonuniform loading and buckled-skin load- 
bearing capability 

Below, we describe a high-accuracy numerical method for analyzing stability of 
built-up composite/metal structures subjected to subcritical loading with compressive and 

shear forces which vary along the cross section contour. 

The structure is a cylindrical or slightly conical system of interacting rods, plates 
and shells made of composites or metals and subjected to mechanical and thermal loads. 

Account is taken of variable curvature, loads, thickness, initial out-of-plane 
displacements, and subcritical deformation of the system along the contour. Also, the skin is 
assumed to be multi-layered; elements, anisotropic; stiffeners, discrete; and global/local 

buckling, likely. 

The method is based on 

formulating the corresponding homogeneous boundary-value problem, 

— separating the problem variables, and 

solving numerically the stability equations for the structural elements by means of 

the discrete orthogonalization method of [10]. 

Of course, all the structural components are linked into a system by using finite element 


methods. 

Consider an isolated anisotropic non-circular cylindrical shell that in the subcritical 
state is loaded with stress resultants N x (s,t) , Nj(s,t) , and N xs (s,t) (Fig. 2,8) where t is the 
load growth parameter, t>0. General equations of structural stability which describe 
buckling of the shell may be represented as three groups of relations (refer to [4]). 

I. Neutral equilibrium equations: 



jQ L _gQs __^x dv Ng fgv | w'j 

r " ax as R ax r las rJ’ (2.20) 
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II. Elasticity relations (from [11, 12]): 

f{N}]_ [[B] [C]lf{e}l 
|{M}J LtC] [D]J\{ae}J ' 

Here, the meaning of new symbols is as follows. 

{N}=|_N x N s S_] T , {M} =[M x M s hJ T , {e} =L e x e sYxsJ > 
{ae} = |_ae x ae s ae xs J T , [B] = [b u ] = [B] T , [C] = [c y ] = [C] T , 
[D] = [ D y ] = [D] T , i,j = l,2,3. 

III. Geometric interrelations: 

du 5v w / o fl V 

S *”as + R Zi' 

v = ^ + * + (e° + e„k, 
y xs as ax v s U - 

ae x _dQj_ _*!_**'» 

se x=^r- *»-& 0S R0X’ 


5w Q 5w v. L 0 o = _^L + x! 

0 x “ ax’ 05 Ss + R’ 1° ds ds R 


'o Q o 


Hereafter, 

u , v, w , 0 X , and 0 S are additional linear displacements and rotation angles of 
the normal to the skin; 

g g s) y xs , ae x , ae s ,and ae xs are strains and curvature variation components, 

w 0 (s) and w°(s) are initial and "accumulated" subcritical out-of-plane 
displacements; and 

By, C y , and D y (i,j=l,2, 3 ) are shell stiffness parameters. 

Boundary conditions for edges {s=0} and {s=b} of the isolated shell are described by, 


w| s= o0-yo) + Qsl s=o T'o=°’ 

w| s=b (l-Yi) + Qs| s=b Yl -°* 


e s| s=0 ( l_5 o) + M s| s=0 5 0 = 0 ’ 

0sUO-»,) + M # | s . b 8 1 «O, 


u| s=0 ( 1 - ( Po) + S| s=0 9o = 0 > 

u ls=bO ~ S| s=b cpi =0, 

(2.23) 

v| s=0 ( 1 -Vo) + N sl s=0 Vo =0 ’ 

v ls=b( 1_V l>l)+ N s| s=b Vl =0. 



Here, the constants y 0 , 8 0 , cp 0 , Vo> Yr> ^i> an( * ta ^ e the values ® or 1 *° r 
cases with displacement control and load control, respectively. Lastly, additional symbols are, 

— 9u o sx 

s = s - | - N s> — + N S — . Q S =Q S + ax • 

n d\ of dv w) 

n s =n s + n; x - + n?(- +i J. 
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The critical value t, of the load parameter t_ is sought to be the minimum value of t 

at which the homogeneous boundary-value problem (2.20) * (2.23) has a non-zero solution 
(subjected to the extra condition of simple support at shell ends {x=0} and {x=a}). 
Underlined terms in the above equations provide options for analyzing general stability of 
the structure in combination with local and temperature-induced buckling. 

Stability equations for a built-up stiffened cylindrical structure composed of 
interacting cylindrical shells, plates and rods may be written in a similar way. For instance, an 
i-th web (see Fig. 2.8) may be analyzed by employing relations similar to (2.20) - (2.23) in the 
local coordinate system (X, , X 2 , X 3 } which are provided in the work [13]; the latter also 
gives compatibility conditions for displacements in the skin and stiffeners. Stability equations 
for plates in the structure may be derived from (2.20) - (2.23) by assuming that 1/R=0. 

All equations are written for each shell, plate and stiffeners. At all joints of shells 
and plates the equations are complemented with the fitting conditions which include 
equilibrium equations and the condition of generalized displacements being identical at each 
joint (refer to [13]). Free and supported longitudinal edges of some elements must meet the 
boundary conditions (2.23). 

The major difference between the new system of homogeneous differential 
equations and the one in [13] is that the former system includes general elasticity relations 
(2.21) which are characteristic of anisotropic composite structures; this circumstance 
notably complicates solution procedures. 

To begin with, we have to consider a particular problem with no subcritical shear 
load and no orthotropy in elements — more specifically, 

N sx=0> B i 3 =B 23 = Ci 3 =C 23 =D i 3 = E > 23 = 0 . The solution for each shell may be written as, 


{P} = {p n )cos^, { F} = {F n | sin-— , 

a 3 . 

where 

{P} =L U0 x Yxs *xs S HQ X J T , {P n } =L u n Qxn.Yn *xsn S n H n Q xn J T , 

{F}=[v W0 S 8 X e s ae x ae s N x N s M x M s Q s J T , 


(Pn } L v n w n e xn e n ®xn 86 n ^xn ^n^ xn M n Q n J 

In addition, n is the number of longitudinal half-waves in the shell, and the functions with 
the subscript n are only depending on the coordinate s. 

Stability equations (2.20) - (2.22) can in this case be transformed into the system of 
8 ordinary differential equations (at any values of n and t as the work [13] states): 


d{u} 

ds 


[C]{u|, 


(2.25) 
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Here, 


{u} =|w n 0 n M n Q n u n v n N n S n J T is the column vector to be sought, 


S„=S„ + N^, N„=N n+ N? 


dv, 


w ^ 


ds R 


nn 


. Qn=Qn- — H n 


(2.26) 


[c] is the matrix whose coefficients are known functions of s. 

Boundary conditions for an isolated element may be written as, 

[r 0 ]{u(0)} = o, [r,]{u(b)} = o, 

where [r 0 ] and [r,] are constant 4 by 8 matrices. 

The problem (2.25) - (2.26) may be efficiently solved by utilizing the discrete 
orthogonalization method (which requires solving some Cauchy problems). For a point 


{s=b} the solution is, 

(u(b)} =[Z(b)]{C(b)} 

and the condition for obtaining the critical value t* is written as follows: 


det[D(t)] = 0, (2 ' 27) 

where [d] = [r, ][Z(b)] and [Z(b)] is the established solution matrix. 

To analyze a built-up cylindrical system, the following procedure is applicable. 
Solutions in the form of (2.24) are written for all elements. Equations (2.25) are obtained for 
each shell and plate while assuming the unit-displacement boundary condition - for 


example, for the following conditions at points {s k - 0} and {s k -b k }: 
= 1 , ej = u k = v k ck _ =w k =e£ = uj = v k | =0. 


n I s k =0 


Here, the superscript k is the identification number of the element. Proceeding in the same 


way, generalized stiffness matrices are established for each element |[K'] j and the entire 

structure ([K]) ; use is made of conventional algorithms of assumed-displacement finite 
element methods. The resulting equations for neutral equilibrium become linear algebraic 


equations (at any n and t): 

[K]{y}=0, <2.28, 

where (y) is the generalized displacement vector for all nodes in a cross section of the 
system. The vector comprises all nodal displacement vectors: 

{y , } = L W ni ni u ni v niJ > i = l,2,...,N 
The equation for determining t* takes the form, 

det[K(t)] = 0. 


(2.29) 
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Initially, we should use (2.29) to find t n at various values of n; thereafter, t. is sought as 
min t n at n = n». Finally, (2.28) is utilized to determine the eigenvector {y.} and the 

n 

structure buckling shape. 

In the general case the subcritical shear load is non-zero and/or elasticity relations 
(2.21) for at least one element have a general form; then the solution to the structural 
stability problem assumes a form differing from (2.24): 

{P} = { P n } cos y + {p n } sin y , {F} = { F n } sin y + {F n } cos y . (2.24') 

Here, variables of the first group ({P n },and {F n }) are the same as in (2.24), 


whereas variables of the second group are 

j^n } = [_ u n ®xn Y n ^xsn Qxn J » 

{Fn } = L V n w n e xn e n ^xn ®n ^xn M xn M n Qn J * 

In this case the structure is assumed to be rather long; i > 0 denotes the buckling 
half-wave for the longitudinal direction (in practical analyses it suffices to adopt 
t = = a/n , n = 1, 2, 3, ...). 

Further, we can substitute (2.24 1 ) into shell stability equations (2.20) - (2.21) to 
separate the variables and obtain 16 simultaneous ordinary differential equations below (in 

which the prime sign symbolizes differentiation with respect to s; 0° =0° + 0 q) : 


W n = -0n + 


n R ’ 


M„ =(Nj + N?0 s o2 )0 n + Q n +2n n H n + 0,°N„ - NXs„ 


®s^sxFn v n "*"^sxM-n w n> 


N r 


Qn =Pn(N° x +N°0 s O2 )w n -p 2 n 0 s °Njv n +-^ L -p n 0 s °S n 4 
+ Pn®s^sYn + Pn^xn “Pn^xs^n + Pn^xs®s u n> 
u n = Yn — Pn v n + ®sPn w n» 

v n _e n _ ^ ® n > 

N; = -f.?,N?0 s °w„ - ^ + (1 ;(n; + Nj)v n + M n S n - 
- n „ N“r „ + N? x( . „ (s„ - 0 S °0„ ) - NV„ u „ , 

= ^n^.x u n ~ M n^xn “l i n^xs(Yn + M-n v n "®s^n w n)’ 
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w' =-e n + 


R ’ 


— ® n » 

M'n =(N° + N?e s 02 )e n + 5n +0s O ^n - - 2p n H n 


-0s°N»H n v n _ N s X H n w n , 


N r 


5; =M 2 n (N° x +N?0 s O2 )w n V n e s °NSv n +^ L + ^ n e s s n 
-|a n 0 sNsY n +Hn^xn +^nNs X 9 n + ^n^sx®s u n* 


U„ = Yn + M-n^n ~ ^Vn^n- 

~ ^ n Q 0 Q 

v " n 'IT -05 Un ’ 


N; =-^SN?e s °w n -% + ^(nJ + N°)v n -n n S n + H n N$Yn - 

~Ns X p n (e n - 9s®n) - N sx p n u n> 

I; = HnN° x U n + M-n^xn + Fn^sx^ n ~ ^n v n + M-n^s w n)’ 


(2.30) 


Additionally, we introduced \i n =nji n . 

The similar way may be used to transform the equations for rods, the boundary 
conditions of form (2.23), and compatibility conditions for longitudinal edges of elements in 


the system. 

Substitution of (2.24 1 ) into elasticity relations (2.21) allows us to express the 
variables e n , ae n , Y n > £ n • » Yn * N xn » M xn > >N xn , M xn , H n (which are present in 

(2.30)) in terms of basic variables 


{u}=[w n 0 n M n Q n u n v n N n S n w n G n M n Q n u n v n N n S n J T . 


(2.31) 


For example, values of { 8 n , se n , Yn } and { s n , £ n , Y n > should be determined from the 


following systems of equations: 


B 2 2 + Ns 

C22 

B 2 3 

£ n 


N n +Nj x p n v n+ N?0 s °0 n 

C22 

D22 

C 23 


> = ■ 

M n 

S n +N? x p n u n -N°(p n v n -0 s V n w n ) 

B32 

C 32 

B33+N? 

.Yn \ 




B 2 l^n 

“ C 2 |Hn 

C23l i n 

u n 

+ 

C 2 l^n 

- D 2 |Pn 

D 23Fn 

w n • 


_B 3 l^n 

-c 3 ,p5 

C 33Fn_ 

A. 


(2.32) 
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b 22 + n" 
C22 
B 32 


C 22 

D 22 

C 3 2 


B 23 

C 23 

B33+N? 


r~ 


se t 


Un 


N n -N^ n v n+ N“0 s \ 


050; 


M, 


S n -N^ n u n -N;(-H n v n+ e s V n w n | 


- 


B 2 !^n 


^23 Hn 

C 2 lHn 

1^21 M-n 

E> 23 Hn ^ 

_B 3 lHn 

C 3 l^n 

C 33 Hn_ . 


u. 


w . 

e. 


(2.33) 


and, similarly, the stress resultants and moments {N xn , M xn , H n } and { N xn , M xn , H n } 
would be found from the equations, 


N xn 


-B U Pn 

C[[Pn 

-C 13 H 

n 


M xn 

= 

_ C[lHn 

D|)Hn 

-D[ 3 p 

n 

w n > 

A . 


L -C 3 iP n 

D3iPn 

-2D 33 |i n J 

l 

V 

A/ 


B)lHn 

QlHn 

C)3Hn 

Up 


“ 

'Mxn 


CilHn 

D|lPn 

D[3Hn 

W, 


> + 

[ H n 


_C3|P n 

D3iHn 

2D 33 p n 

IV 
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B| 2 

C,2 

C 32 


C, 2 

D, 2 

D 32 
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C 33 
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IV 


(2.34) 


Bj 2 
C,2 

^32 


C,2 

D|2 

D 32 


B, 3 

C, 3 

C 33 


Un 


(2.35) 


The relevant matrix-based problem statement and the solution method completely 
coincide with those proposed above; however, dimensions of all vectors and matrices are 
twice as much, so in the general case 


{yi } =L w "' e ni u ni v ni W nL 0 ni U ni v ni J , 1 - 1 . 2 .....N. 

Integral stiffnesses B ( , C,, and D, of multilayered composite element should be 

evaluated on the basis of elastic properties of layers present in the stack; use is made of the 
methods exposed in [1 1, 12); if early buckling takes place, these are replaced with tangential 
characteristics considered in Section 2.1. Any component in the system can be isotropic, 

metallic or orthotropic "by design" (refer to [1 0]). 

This means that the method for numerically analyzing built-up cylindrical 

structures that had been developed by the authors in [13] is now made capable to deal with 

structures which incorporate composite materials. 

The present method and the respective software programs can be employed to solve 
many classical problems (for which the researchers have provided exact and/or reliable 
approximate solutions), new problems of structural buckling, as well as complicated problems 
on stability of built-up composite structures that appear in designer's practical work. 

For example, Table I compares the new numerically obtained results with the exact 
solution of [12] for a rectangular orthotropic plate simply supported at all four edges. 
Analyses were performed assuming the following parameters: 
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_ the total number of orthogonalization points for each element, n 0 = 10 , 

— the total number of Runge-Kutta steps between orthogonalization points, "p-K = 10 ■ 

— the allowable relative error in the t„ determination procedure, e = 0.001 . 


Table 1 


N x b 2 ^ 
k = Nx - = U 2 s/D,Dj. 


e = 10 3 
n 0 =10 


n P-K 


= 10 



It is clearly seen that the error of the numerical solution at the aforementioned 
procedural parameters does not exceed 0.05% and the critical load thus obtained ts above 

the exact value. 

A similar result is seen in comparison of the exact value of stability coefficient, 
k=7.69 (refer to [14]) and the numerically obtained one, k=7.6916, for the problem of 

buckling of a square isotropic plate with clamped longitudinal edges. 

As for the boundary condition version mentioned, we find it interesting to estimate 
the accuracy of the work [15] simplified solution which is widely used in analysis of 
composite plates. The solution is compared with the newly obtained data in Table 2 an 
shown to overestimate the critical load by 1 .6% to 4/o. 

Table 2 
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Further, we solved the problem on stability of a longitudinally compressed 
composite pane, whose stiffeners are longitudinal webs with a rectangular cross se«,om 
Panel dimensions are a=500 mm and b=300 mm; stiffener spacmg ,s 60 mm, s 
thickness x height, 3.64 mm x 30 mm; and skin thickness, 2.34 mm. s 1 en 
stacking sequence and the skin, P* charactensttcs 

E, =20930 kg/sq.mm, E 2 =800 kg/sq.mm, G, ,=520 kg/sq.mm, v„=0.26, 5=0.13mm. 
Consideration was given to 

_ global buckling (Fig. 2.9a) for the panel with unsupported longitudmal edges 

a " d - local buckling (Fig. 2.9b) for the panel with simply supported longitudinal edges. 

In the case of global buckling the panel has one half-wave ( n . = 1 ) and the cnt.cal 
compressive load P. = 57500kg; and in the case of local buckling, n.=5, 
l. = a /n. =100 mm, and P. =92500 kg. 

Also, we treated the problem on buckling of a shear-loaded long composite plate 
(a=2000 mm’, b=500 mm) whose layers have E, =27000 kg/sq.mm, E,=2000 kg/sq.mm, 
<3,2=980 kg/sq.mm, v 12 =0.22, 6 ply =0.1 mm. The stacking sequence was defined as 

[0/90/ct4]s where the angle a was varied from 0 to +90 and from 

Figure 2. 10 represents the dependence of the critical shear load N s< on the angle o 
of the eight middle layers. One can see that if the angle a is negative, the critical shear load 
is maximum at the angle of approximately -60°. If the angle a is positive, the cnt.cal shear 

load is much lower. 

Figure 2.1 1 demonstrates how the critical compressive load N, of the same p ate 

depends on the angle a of the eight middle layers. The critical load is maximum at a =45". 

Of particular interest is stability of a composite plate in which certam layers are 
disbonded over a domain with the length ? . Even if the initial layup was symmetnc, the 
debonding makes the panel become a structure with two asymmetric stacks over the dotnatn 
mentioned. This system may well be analyzed by the above methods and software programs 
since these are capable to treat the structure as a system of asymmetric elements included m 

the general model. 

Used as an example may be the problem for a composite panel with the stacking 
sequence [90/0/+452/-452/90,/0],. Each ply is assumed to have the following characterist.es; 
£, = 18000 kg/sq.mm, E 2 =900 kg/sq.mm, G 12 =514kg/sq.mm, v, 2 =0.31 , S P , y =0.05 mm. Panel 
sizes are a=120 mm and b=40 mm. Now we presume that three upper layers at the center of 
the plate are not bonded over the length 4 which is varied in the course of the analyses. 
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Theoretical results are shown in Fig. 2.12 together with the analytical model 
comprising five elements and five nodes. Elements 2, 3, and 4 feature an asymmetric ayup 

and anisotropy. 

In the debonding domain the subcritical stress resultant is distnbuted between 
upper and lower stacks in accordance with the deformation compatibility condition for t e 


limiting points of the debonding domain. 

Dependence of the critical load N, on the relative debonding domain length 4 is 
shown in Fig. 2. 1 2 by the solid line. The dashed line represents the critical load for an intact 
plate, and the dash-and-dot line corresponds to the critical load for a stepw>se plate wit 

disbonded layers removed in the domain with the length t, . 

Analyzing the solid line in Fig. 2.12, it is clear that a short debonding domain 
(^/b<0.1) does not cause the substacks to separate, so the plate is buckling as a whole 
following the usual global buckling shape. One can see an insignificant decrease in the 
critical load for the interval {0<^/b<0.1}; the reason is the stack stiffness being reduced m 


the debonding domain. 

If %lb>0.\ the thin part of the plate is buckling, with the rest of the plate being 

almost straight and capable of carrying a load increment. At higher loads the stiffness of the 
thin part gets reduced, which is followed by global buckling; the critical load at the moment 
of reaching the load-carrying capacity limit corresponds to the critical load for a panel with 
a weakened part over the disbonded area. 
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CONCLUSION 

In compliance with the NCCW- 1-233 the investigation into the influence of the skin 
postbuckling load-bearing capability on general stress and strains and on stability of 

composite structures has been accomplished. 

The following results have been generated: 

_ a simplified method for analyzing stresses and strains and load redistribution 

between structural components which is caused by skin buckling; 

_ a nonlinear finite element method taking into account the skin postbuckling 

load-bearing capability; 

_ an analytical solution to stability problems; evaluation of the effect of the skin 
postbuckling behavior on general stability of stiffened composite panels and shells; 

— a high-accuracy numerical method for analyzing stability of built-up cylindrical 
structures that allows for discrete stiffening, nonuniform multiaxial load, and skin 
postbuckling behavior. 

The results and proposals on further research work would be discussed with NASA 


experts. 
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Figure 1.8. Maximum stresses in lower panel skin: load B 
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Figure 1.11. Mean stresses in compression-loaded panel skin and 

stiffeners: load A 



Figure 1.12. Mean stresses in compression-loaded panel skin and 

stiffeners: load B 






Figure 1.14. Maximum stresses in compression-loaded panel skin: load B 
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Figure 1.15. 
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Figure 1.18. Normal stress diagram 
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Figure 1.19. Normal stress in stiffeners at clamped end 
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Figure 2.8. Element for evaluating anisotropic cylindrical shell 
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Figure 2.9. General and local panel buckling 
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Figure 2.12. Stability and load-carrying capability of anisotropic 
composite plate with partial debonding 




